МІНІСТЕРСТВО ОСВІТИ І НАУКИ УКРАЇНИ
Національний університет “Львівська політехніка”
Інститут післядипломної освіти
ЗВІТ
Про виконання лабораторної роботи №2
«Чисельні методи розв’язування системи лінійних рівнянь»
з дисципліни «Чисельні методи»
Тема роботи: Чисельні методи розв’язування системи лінійних рівнянь.
Мета роботи: Ознайомлення на практиці з прямими методами розв’язування систем лінійних алгебраїчних рівнянь.
1. Теоретичні відомості
Система лінійних алгебраїчних рівнянь
Математичні моделі багатьох технічних задач представлені системами лінійних рівнянь. Багато методів розв’язання нелінійних задач також зводяться до розв’язання деякої послідовності систем лінійних алгебраїчних рівнянь (СЛАР).
Прямі (точні) методи дозволяють розв’язати систему рівнянь за скінчене число арифметичних операцій. Якщо всі операції виконуються точно (без помилок округлення), то розв’язок заданої системи також отримуємо точним. До прямих методів належать: метод послідовного виключення невідомих (метод Гауса та його модифікації: метод головного елемента, метод квадратного кореня, метод відображень та ін.), метод ортогоналізації, метод LU-розкладу.
Системою лінійних алгебраїчних рівнянь (СЛАР) називають систему вигляду
,
(1.1)
де – невідомі, – вільні члени системи, – коефіцієнти системи.
У матричному вигляді система прийме вигляд:
(1.2)
де , , .
(1.3)
Розв’язком системи (2) називається n-компонентний вектор-стовбець
, який перетворює матричне рівняння (2) у вірну числову тотожність.
Система може мати один розв’язок, безліч або жодного. Кількість невідомих m в системі називають порядком СЛАР. Систему лінійних алгебраїчних рівнянь називають:
- сумісною, якщо вона має хоча б один ненульовий розв’язок;
- несумісною, якщо СЛАР не має жодного розв’язку;
- визначеною, якщо вона має тільки один розв’язок (випадок, коли m=n);
- невизначеною, якщо система має безліч розв’язків (m(n);
- виродженою, якщо головний визначник системи дорівнює нулю;
- невиродженою, якщо головний визначник системи не дорівнює нулю.
Дві системи називаються еквівалентними, якщо ці системи сумісні, визначені і мають однаковий розв’язок.
Систему лінійних алгебраїчних рівнянь можна розв'язати на ЕОМ числовими методами, якщо вона сумісна, визначена, невироджена.
Блок-схеми алгоритму Гауса з вибором головного елемента (рис. 1.1) та методу LU-розкладу (рис. 1.2):
Рис. 1.1. Схема алгоритму розв’язання СЛАР методом Гауса з вибором головного елемента
Рис. 1.2. Схема алгоритму розв’язання СЛАР методом LU-розкладу
2. Хід роботи
Завдання (варіант 5). Написати програму розв’язку системи лінійних алгебраїчних рівнянь у відповідності до варіанту:
1) методом Гауса з вибором головного елемента;
2) методом LU-розкладу.
Вимоги до звіту
Звіт до лабораторної роботи повинен містити такі структурні елементи:
Титульний аркуш.
Тема.
Мета.
Короткі теоретичні відомості.
Алгоритм розв’язку СЛАР.
Текст програми з коментарями.
Вигляд реалізованої програми.
Висновки.
Вимоги до програми
Програма має передбачати наступні можливості:
Автоматичний розв’язок СЛАР відповідними методами.
Ввід вхідних даних вручну: задати порядок СЛАР, елементи матриці коефіцієнтів та вільні члени системи.
Передбачити можливість некоректного введення даних.
Передбачити вивід покрокового виконання для кожного методу.
3. Текст програми методу Гауса з вибором головного елемента
#include <stdio.h>
#include <conio.h>
#include <clocale>
#include <math.h>
int const n = 3;
//якщо 0 - то ввід вхідних даних з клавуатури
int const testMode = 1;
void mainElem( int k, double mas[] [n + 1], int n, int xx[] ) {
int i, j, i_max = k, j_max = k;
double temp;
//Шукаємо максимальний по модулю елемент
for ( i = k; i < n; i++ )
for ( j = k; ...
